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Abstract 


In  this  report  we  explore  the  problem  of  obtaining 
confidence  intervals  and  tests  of  hypotheses  for  the 
parameters  of  a  new  family  of  life  distributions  derived 
In  [1].  Several  statistics  whose  distributions  are 
Independent  of  the  shape  parameter  are  investigated  in 
terms  of  accuracy  of  confidence  bounds  for  the  distribution 
scale  parameter.  Percentage  points  of  the  statistic  yielding 
the  most  accurate  confidence  bounds  on  the  scale  parameter, 
among  those  investigated,  are  tabulated  for  samples  of 
size  n,  n*2(l)5. 
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0.  lilt  roilurt  Ion 

In  (l)i  Birnhutim  and  S;mnik,r'>  derived  a  new  family  of  Ilf*' 
distributions  from  plausible  nmt.  IdiT.it  lout*  of  the  physical  behavior 
of  fatigue-  crack  grnwt li  undt-r  repeated  loading,  Certain  point 
estimators  for  tin*  par a  me  tors  of  this  family  of  d  1st rl but  ions  v«- re 
investigated  in  (2).  In  t lie  following,  wo  explore  the  problem  of 
obtaining  confidence  intervals  nnd  tests  of  hypotheses  for  these 
parameters. 

We  denote  by  £F(u,f. )  tne  two-parameter  distribution  of  a  non- 
negative  random  variable  derived  in  (I)  And  defined  by 

(0.1)  W[-  Ut/fc>)  for  t  >  0 

a 

where  a  >  0,  8  >  0  and 

(0.2)  C(t)  -  St  -  ±  , 

/t 

and  9?  is  the  distribution  function  of  the  standard  normal  variate. 

The  parameter  a  determines  the  shape  of  the  distribution  while  the 
median  8  is  a  scale  parameter.  Iterative  numerical  procedure#  for 
the  computation  of  the  maximum  likelihood  estimates  of  a  and  r  are 
developed  and  a  simple  estimator  c  of  8  is  derived  in  (2).  We  show 
here  that  the  maximum  likelihood  estimators,  and  the  estimator  r  as 
well,  even  though  they  have  optimal  asymptotic  properties,  are  such 
that  they  cannot  form  the  basis  for  invariant  tests  and  confidence 
nleivuls  for  either  of  the  parameters  for  all  sample  sizes.  Several 
st  at  I  *•  t  I  <  s  whose  distributions  are  independent  of  the  shape  parameter 
u  are  investigated  In  terms  of  their  consequent  confidence  regions 


Im  ill*-  in  .1 1  *•  |mi  .ii'k-i *m  i‘ .  Tin*  percentage  point#  of  the  moat 
.»*  <  i'i't .it- !*•  *>t  .it  lot  ii’»  among  tlmae  Invent  1 gated  arc  evaluated  by 
Mi'iii «  Carlo  nu*tlii*i|n  ami  the  remitta  presented.  Further,  some 
i  ompnr  Umit,  of  power  are  alao  obtained  by  the  same  tncana  and  the 
curve*  exit ll'i ted. 
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1 .  Maximum  Likelihood  Estimators 

In  order  to  define  the  maximum  likelihood  estimators  of  a  and  p 
we  first  define  for  a  given  set  of  positive  random  variables  T^,...,Tn 
the  arithmetic  and  harmonic  means,  respectively 


and  the  harmonic  mean  function 


Then  if  . Tr  represents  a  sample  of  independent  random  variables 

each  distributed  as  ^(a,8),  the  maximum  likelihood  estimator  8  of  8 
is  the  unique  positive  solution  of  the  random  equation  g(x)  *  0  where 


g(x)  -  x2  -  x[2R  +  K(x)]  +  R[S  +  K (x) ] . 


The  maximum  likelihood  estimator  a  of  a  is  then  given  in  terms  of 
8,  S,  R  and 


These  results  were  given  in  ( 2 J . 

The  simplified  estimator  8  of  8,  which  is  shown  in  [2]  to  be 
consistent  and  for  small  values  of  a  virtually  the  same  as  the  maximum 
likelihood  estimator  8  is  given  by 

8  ■  /SR. 
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If  tin*  sample  si/.f  is  n  =  2,  Llu*n  it  nl.su  follows  from  tin* 


results  of  (2]  that  tin*  maximum  likelihood  estimators  of  a  and  R 
sat  i  s f  y 


Wo  note  that  for  n  =  2,  the  maximal  invariant  for  confidence  sets 
and  tests  of  hypotheses  independent  of  a  is  the  r  it  io  (T.,/.- )/r  (T  . /,•  ) . 
For  n  =  2  there  is  apparently  no  function  of  B / R  and  a  to  which 
this  ratio  is  equal.  One  would  infer  that  functions  of  a  and  R  (or 
cannot  he  used  to  obtain  confidence  sets  .for  R.  In  fact,  we  now  demon¬ 
strate  that  for  n  =  2  one  cannot  have  confidence  sets  and  tests  of 
hypotheses  for  either  a  or  R  based  on  functions  of  a  and  B. 

Let  denote  a  random  variable  with  law  f7(u,l).  In  distribution 

we  see  that 


(1.1) 


1 


since  each  is  equal  in  distribution  to  its  own  reciprocal.  This 

property  we  call  the  reciprocal  property.  Thus  (1.1)  says  that  R./i 
has  the  reciprocal  property.  We  now  find  the  distribution  of  the  product 
XjX2 ,  call  it  G: 


G (t )  =  P[X1X2  <_  t]  -  P[XX  1  tX2] 
00 

=  /  c(tx)]d??[~  £(*)]• 

Jq  a  a 

Letting  p  =  f,  ^  we  see 


Thus  in  terms  of  the  distribution  of  G,  which  contains  a  as 
a  nuisance  parameter,  we  can  express  via  (1.4)  and  (1.3)  the  distributions 
of  a  and  g.  Note  how  the  parameter  a  is  "scrambled"  in  the  distribution 
of  a,  in  fact,  it  appears  to  be  almost  inextricably  bound. 

Using  the  fact  that 

£(t)  =  2  sinh  (?-~2  " )  for  t  >  0, 

an  identity  given  in  [3],  and  some  straightforward,  but  tedious,  algebra 
one  finds 

(1.5)  ~  Ut^(ay)]  =  ~[slnh  ~  +  sinh  1  (^)] 

=  y  sinh  (■!=■)  +  /(- ) 2  +  y2  cosh  (-|) . 
z  a  z 

Thus  we  <  on elude  that  for  exact  inference  the  distribution  of  3,  based 
on  two  observations,  would  be  of  use  only  if  a  were  known.  Since 

“■  UttKay)]  1  y  sinh  (-j)  +  { y  [  cosh(y)  =  2(t,y) 

say,  we  have 
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(1.6) 


G(t)  i  f  9lU(t,y)]d9f(y). 

—  00 

Thus  finding  a  tj.  such  that  the  right-hand  side  of  (1.5) 
equals  l-e  would  enable  one  to  obtain  a  conservative  lower  confidence 
bound  for  8  (which  is  the  one  of  primary  interest  in  life  studies) 
since  by  (1. 3) 

(1.7)  p[(e/^7>  13]  -  G(t1-e)  >  l-e. 

This  bound  is  of  level  l-e  and  invariant  with  respect  to  a.  Whether  or 
not  it  is  of  use  would  depend  upon  the  outcome  of  numerical  tabulation  and 
the  practical  situation. 

The  conjecture  that  for  all  samples  of  size  greater  than  two, 
the  maximum  likelihood  estimators  of  both  a  and  §  would  contain  the 
nuisance  parameter  a  would  seem  to  be  substantiated.  At  the  same  time 
a  result  as  explicit  as  (1.2)  and  (1.6)  appear  to  be  a  formidable  task 
for  large  sample  sizes. 
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2.  Invariant  Confidence  Intervals  and  Tests  for  8 


Suppose  we  have  a  sample  T^,...,Tn  ot  independent  £7(ci,B) 
variates.  From  the  equations  (0.1)  and  (0.2)  it  follows  that 

T 

Zi=^-C(y-)  for  1-1,... ,n 


are  independent  standard  normal  variates.  Maximal  invariants  for 
estimators,  tests  and  confidence  bounds  for  6,  invariant  with  respect 
to  a,  are  functions  of  ratios  of  such  variates.  We  can  try  to 
utilize  well-known  results  for  the  normal  distribution  to  obtain  bounds 
on  B  functionally  independent  of  a . 

First  let  us  consider,  for  k  +  m  “  n,  the  statistic 


Z 


2 

i 


which  has  Snedecor's  F-distribution  with  k  and  m  degrees  of  freedom 
in  numerator  and  denominator,  respectively. 

From  standard  tabulations  we  can  find  6^,  &2  such  that 


P[62  ‘  Uk,n  ‘  Sl>  ■  Y 


for  any  prescribed  ye  (0,1). 

The  question  is  whether  these  percentage  points  of  the  F-distribution 


can  be  used  to  obtain  confidence  bounds  on  the  unknown  parameter  6 

2  1 

independent  of  a.  Using  the  fact  that  £  (x)  =  x  +  —  -  2  we  study  an 
observed  value  u,  (B )  as  a  function  of  B  and  obtain 

fC  jlTl 


U.  m(S) 
K  ,m 


1  k  7 

k  £  5  (tl/6) 

K  i=l 


sx  +  (B  /r x )  -  26 


i  i  s2  +  - 26 

m  k+1 


i  k  / 1  k  -A-1 

where  s  =  r  ^  t.,  r.  a  r  ^  t.  )  are  the  nr i  Mimetic  and 

1  k  T  1  1  \k  w  1  / 

harmonic  means  of  the  first  k  observations,  with  s^,  r ^  defined 


similarly.  Consider  the  quadratic  equation  representing  either  the 


numerator  or  denominator.  One  sees  that  it  has  roots  r,  l  1  t  /I  -  — 

\  ri  / 

but  from  the  well-known  inequality  between  r,  s  we  infer  that  both 
roots  are  complex.  Therefore,  no  real  roots  or  poles  exist  for  the 
rational  function  g(8)  for  8  >  0.  Unfortunately,  closer  examination 
shows  that  the  distribution  of  m  is  virtually  invariant  with 
respect  to  8  so  that  inversion  to  obtain  meaningful  confidence  bounds 
is  not  possible.  To  see  this  note  that 


uk,m(0)  ’  (S1/S2)>  uk,.W  '  (r2/rl) 


which  both  have  an  expected  value  near  unity.  Now 


sgn  u^  ^(6)  »  sgn  (aB  +  b8  +  c) 


11  S2  S1 

where  a  - - ,  b  - - ,  c  ■  s.  -  s_.  Thus  the  function 

r2  ri  rl  r2  12 

u,  (8)  has  one  local  maximum  and  one  local  minimum  but  the  actual 

k,m 

behavior  depends  upon  the  values  of  a,  b,  c.  Note  that  a  will  always 
be  near  zero  so  that  the  quadratic  equation  is  nearly  linear.  The  roots 


b  i  -4ac  _  |  b  | 


7  c  r  ac 


(sgn  b  ±  1)  +  fTT  + 


+  0(a) . 


Examination  shows  that  one  of  the  roots  is  near  |^,|  for  |  a|  small. 
Some  actual  computer  runs  on  the  IBM-360  were  made  for  the  case 
k  =  m  =  50  and  yielded  the  following  values: 
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SUp  Uk,m(B) 

inf  u 

k,i 

1.1941 

.5968 

1.8273 

.8473 

1.3341 

.9117 

The  range  of  values  was  very  small.  This  idee  was  then  abandoned. 
Consider  the  t-distribution  with  n  degrees  of  freedom 


V 

n 


i=l  /n 


Again  using  the  definition  of  the  function  £  we  examine  an  observed 
value  as  a  function  of  B,  namely 

££( t./B)  *//np  -  B*^ 

vn(e)  =  1 . . . 2-9 — T 

/zc2( ti/B)  (s  +  -  23)^ 

where  s,p,q,r  are  the  means  of  (t^,...,tn)  of  order  1,  ,  -'5,  -1 , 

respectively. 

One  checks  that 


sgn  v'(B)  =  sgn [&(•'  —  -  ^2-)  +  /npT  -  st-] 
n  q  r  q 

and  thus  v’  changes  sign  only  once.  For  certain  cases  v  (3)  is  monotone 

n  n 

in  B  but  again  the  total  va'iatior.  is  very  small.  Consequently,  we 
see  that  invariance  with  respect  to  a  virtually  forces  an  invariance 
with  respect  to  B,  at  least  in  the  cases  examined.  We  now  try  to 
allow  some  dependence  of  our  statistic  upon  a  for  finite  sample  sizes. 
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Hence  we  seek  a  statistic  whose  distribution  is  invariant  with 


respect  to  B,  depends  upon  a  for  finite  sample  sizes  but  does  not 
depend  upon  a  asymptotically.  Let  us  try 

I  s(t J&)/[  I  ^(Ve)]54 

i-1  1  i-1 

where  $  =  /SR  is  the  mean  mean,  i.e.  the  geometric  mean  of  the 
harmonic  and  arithmetic  means  Our  rationale  for  this  choice  is 
that 

(a)2  -if  £2(T  /8) 

n  i-l  1 

2 

is  the  maximum  likelihood  estimator  of  a  if  3  is  the  maximum  likeli¬ 
hood  estimator  of  6.  For  smell  a  we  know  by  previous  work  that  the  mean 
mean  is  the  Same  as  the  maximum  likelihood  estimator  but  much  easier  to 
compute.  This  we  put  in  the  denominator  and  obtain 


(a) 


2 


c2(t  ±/B) 


Since 


I 


i.j-i 


+  0(i> 


we  see  that  as  n  -+•  “  with  probability  one 

«<4>  -  cd +  v  “ 0,2 2°2  • 

Thus  the  proposed  statistic  cannot  be  invariant  with  respect  to  o  even 
asymptotically.  Consequently,  we  propose  a  modified  statistic  which  has 
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a  distribution  independent  of  £  and  asymptotically  independent  of  a, 
namely 


n 


The  coefficient  has  been  chosen  to  make  S/R  -  1  an  unbiased  estimator 
2 

of  a  .  Consider  an  observed  value  as  a  function  of  £ 


It  is  clear  that  as  a  function  of  B  it  is  not  only  monotone 
decreasing  but  has  a  range  of  «  to  Thus  it  could  be  inverted 

to  obtain  confidence  intervals  for  £  if  percentage  points  were  tabulated. 
Unfortunately,  simulation  by  Monte  Carlo  procedure  of  the  distribution  of 
Yn  revealed  that  although  the  statistic  is  indeed  asymptotically  inde¬ 
pendent  of  a.  the  independence  is  not  exhibited  for  samples  of  size  50 
or  smaller. 

Z2 

We  now  consider  the  ratio,  letting  k+m=n ,  of  W  =  where 
n  T  k  T  1 

Z  *  —  f  £  (•—) ,  Z  *  —  S  £  (~) ,  and  W  has  the  distribution 
1  &  i-k+1  3  2  A  i-1  6 

of  the  ratio  of  two  independent  standard  normal  variates.  This  distribution 
is  the  Cauchy  distribution  F  given  by 


F(x)  *=  4  +  -  arc  tan  x  -«  <  x  < 


Alternatively,  letting  F.  .  denote  Snedecor's  F-distribution  with  one- 

1  |  I 

degree  of  freedom  in  numerator  and  denominator,  we  have 
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7  9  2 

P[-u  •  W  <  u]  =  P[W  <  u~]  =  F  ,(0 

=  F (u)  -  F'-u)  =  2F(u)  -  1. 

The  last  equality  follows  since  -W  has  the  same  distribution  as  W 
Thus  tables  of  the  F  distribution  could  be  used  as  percentage 

point s . 

Consider  an  observed  value  of  W  as  a  function  of  3s 


AT 


/m  r 


w(6)  =  (b/~  -  /kp„)/(3f  --n 

q2  z  ql 


mp  ')  . 


We  notice  that 


w(0) 


fk^2  ,  s  m 

'  m  p  ’  W  °°  ‘  m  q. 


and 


sgn  w(B)  =  fmk  sgn( 


ql  q2 


Thus  the  function  is  either  everywhere  increasing  or  everywhere  decreasing, 
and  we  can  distinguish  exactlv  two  cases: 

Case  1.  pxqL  >  P2q2 
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<'ase  I  !  . 


Pi ‘I  I  P2q2 


wi/il 


Figure  2. 

Suppose  it  is  known  that  for  some  c  we  hav^  P[t  <  W  <  u]  ■  1— c ; 
then  {.•  >  0;i  w(S)  <  u}  is  a  1-e  level  confidence  region.  By 
inspection  we  see  t hat  any  set  of  the  form  {3  >  0:i  <  w(B)  <  u}  could 
consist  of  two  critical  regions,  an  upper  critical  region  or  a  lower 
critical  region  only,  or  the  empty  set  which  would  occur  at  random 
depending  upon  t he  values  of  re^®^^,on  to  l,  u. 

This  is  an  unsatisfactory  state  of  affairs  since  we  should  like 
not  only  to  he  able  to  determine  one-sided  confidence  intervals  but  to  be 
able  to  place  half  of  the  confidence  that  6  lies  in  the  upper  critical 
region  and  half  that  it  lies  in  the  lower.  Moreover,  as  a  matter  of 
practical  convenience,  we  should  prefer  that  one  calculation  suffice  to 
determine  both  upper  and  lower  confidence  intervals.  To  this  end,  we 
propose  to  find  a  method  of  separating  these  regions.  One  realizes 

that  this  is  equivalent  to  ordering  the  observations. 
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3.  Tests  and  Confidence  Sets  for  3  Based  on  Ordered  Observations 


We  observe  that  ^/3  <  1  implies 


x,-  =  UT,  „/3)  =  (T  /3)'*  -  (B/T.  V  <  0 


l  ,n 


l  ,n 


i  »n 


i,n 


and  T.  /g  >  1  implies  X.  >  0.  We  observe  also  that  £(T.  /0) 

i  i  ,n  l  ,n 

is  an  order-preserving  transformation.  We  note  that  for  a  sample  of 

size  2,  Prob[X  <  0]  =  ProbfX,  >  0]  =  (1/2)2.  Hence  Prob[0  >  T„  ]  = 
^•»n  1  ,n  z  ,n 

2 

Prob[$  <  T  ]  =  (1/2)  .  We  can  therefore  determine  iteratively,  by  Monte 
l  fn 

Carlo  procedures,  a  value  (c )  such  that  the  joint  event  X«  -/X  „  >  u(t) 

Z %Z  19Z 

and  X^  2  <  0  occurs  with  probability  e.  Then,  an  upper  confidence  bound 
on  0  at  level  l-e  is  given  by  max[T2  2>g],  where  g  is  equal  to 

(u(c)T^  2-T2  2^^u^e^T1  2~^T2  2^'  Because  of  the  symmetry  of  the  normal 
distribution,  a  lower  confidence  bound  on  0  at  level  l-e  is  given  by 

min[Ti,2,g' [u(e)  ,Ti,2,T2,2]  1  with  g'  =  [u(e)T2  2]/[u(e)/T2  ^l/^  2^* 

Clearly,  because  Prob[X.  >  0]  *  (1/2 )n ,  procedures  based  on  the 

i  ,n 

assumption  X.  >  0  or  X  <0  will  encounter  difficulties  if  n  is 
itn  n,n 

large,  unless  l-e  is  extremely  close  to  1.  One  might  also  consider  a 
test  based  on  the  joint  event  X2  2/X^  2  >  u(e)  and  2  >  0  or  the 
joint  event  X^^  2/X2  2  >  u(c),  X2  2  <  0. 

For  samples  of  size  n  one  can,  therefore,  determine  l  (e )  such  that 
the  probability  that  the  event 


nS(0)  =  V  X  <  0 

i=l  1,n 


and 


W,  (B)  -  V  X.  /X.  >  Kc) 

1  \^i)A  *•"  1-n 


occurs  is  e.  Then  an  upper  confidence  bound  on  0  at  level  l-e  is  given 
by 
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n  ■  max 
n 


2  T*f  /  V  (l/T,  V1,  hU(c),T.  . T  J 

A  *»n  j-i  J»n  1*n  n»n 


(3.1) 


where  hU(e),T.  ,...,T  ]  -  (i/n^T  Kc)!*?  -  2  T)  J/ 

l,n  n  n,n  l,n  1»n 

[/n-1  1(c)/ t!?  -  V  (1/T.  J*1).  From  considerations  of  symmetry, 

l*n  £ 2  1*n 

one  obtains  a  lower  confidence  bound  on  3  at  level  1-e  as 


u  -  min(  2  J  2  (1/Ti  _ )  i  h'U(c),T  ] 

n  i-1  1»n  J-i  J'n  1*n  n*n 


(3.2) 


n-1 


.]/ 


where  h'U(e).T  . T  ]  -  [^=T  £<e)tJ  -  2  T*  J 

1»R.1  n*n  n*n  i-i  i'n 

[i^ri(t)/xjn  •  2  u/t4  nn. 

n,n  jti  J»n 

A  test  of  the  hypothesis  H^:  S  >_  Bq  versus  H^:  B  <  B  based 
on  W^(0q)  rejects  at  significance  level  c  if 


n  n 

5(T1/B0)/C(T1/B0)  >  1(e)  and  nS(8Q)  -  j  ^/Bq)  <  0; 


a  test  of  H2:  B  <  Bq  versus  :  6  >_  Bq  based  on  Wn: 


W  -  "l  c(T1/e0)/?(tII/60) 


rejects  H2  at  significance  level  c  if 


Wn(B0)  ii(c)  and  nS(BQ)  >  0. 


A  two-sided  test  at  significance  level  2e  rejects 
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1,3:  ^  i  ^2  versus 


iiS  . )  >  0  and  W  (3 

1  n 

In  like  manner  one 
confidence  level  l-2c 


H.  :  3  <  3.  or  3  > 
4  1 

)  >  4(c)  or  nS($1) 

can  obtain  two-sided 
as 


<  0  and  W^(3^)  >  fc(e). 

confidence  bounds  at 


w  <  3  <  n  • 

n  n 

Ordinarily,  however,  since  3  is  the  median  of  the  distribution,  one 
would  want  only  a  lower  bound  on  3  or  would  want  to  test  only 
versus 

Clearly,  there  are  many  other  invariant  functions  of  the  ordered 

X^s  that  one  might  consider.  It  is  felt,  however,  that  for  very  small 

sample  sizes  and  relatively  small  values  of  e,  it  is  appropriate  to 

base  tests  and  confidence  sets  on  W,  and  W  ,  with  constraints  either 

i  n 

on  nS  or  on  X,  and  X 

l,n  n,n 

In  the  following  section  a  Monte  Carlo  Investigation  is  made  of  the 

power  of  three  tests  based  on  W,  and  W  . 

1  n 
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4.  Study  of  Power  of  Three  Proposed  Invariant  Tests  Concerning  8 

The  Monte  Carlo  study  described  below  was  designed  to  determine 
which  of  three  proposed  procedures  would  provide  most  powerful  tests, 
based  on  very  small  complete  samples,  concerning  the  parameter  8- 
All  three  procedures  depend  upon  and  W^,  but  these  statistics 

are  combined  differently  in  each  case  with  constraints  on  X^  n 

n 


or  nS  =  £  X  . 

i=l  1,n 

Two  of  the  tests  are  generalizations  of  the  first  two  tests  in 

Section  3  applying  to  a  sample  of  size  2.  The  third  is  the  test 

described  in  Section  3  based  on  W  or  W.  and  constraints  on  nS. 

n  l 

The  first  Monte  Carlo  simulation  procedure  described  below  is 
designed  to  determine  critical  values  for  testing  either  8  ^  8q 

versus  H2:  8  <  8Q,  H2  versus  1^,  or  Hj!  Bj  1  6  i  Bj  versus 
H^:  8  <  8^  or  6  >  82  at  sign^icance  level  e.  In  the  study  we 
considered  samples  of  size  n,  n=2,3,4  and  5  and  values  of  t 
equal  to  .005,  .01,  .025,  .05,  .10.  The  number  N  of  Monte  Carlo 
samples  generated  was  10,000. 

For  each  n,  N  samples  X^,...,Xn  of  independent  stsndard  normal 
variates  were  generated  using  the  method  of  Marsaglia  and  Bray  [A],  For 
each  of  the  N  samples,  the  values 


nS  «  X,  +•••+  X  ,  X  and  X 

1  n  l,n  n,n 


were  determined.  These  values  were  then  used  to  calculate 


nS-X 


W 


nS-X 


1  ^X 


and  W 


n,n 


l,n 


n  ^  x 


n,n 
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The  N  values  of  X  ,  W,  and  W  allowed  us  to  determine 

n,n*  1  n 


estimates  of 


D^z) 

D2(z) 

D3(z) 


;  <  o, 

W.  > 

z] 

n,n 

1 

<  0, 

w  > 

z] 

n,n 

n 

P[nS  <  0,  W2  >  z] 


so  as  to  produce  a  table  of  5^(0  such  that  D^U^Ce)]  ■  e.  The  entries 
of  the  tables  were  obtained  by  letting  (j ) , . . .  »Xn(J )  be  the  J**1  of  the 
N  samples  of  size  n  and  calculating,  e.g.  with  k  ■  1, 

V*>  -i?  £  !Xn,nU>  <  Wl«>  >  *}• 

Here  { A}  is  the  indicator  of  the  event  A,  equal  to  one  if  true 
and  zero  otherwise.  The  valuas  D^(Zj),  were  computed 

from  independent  samples  vhich  form  a  sequence  of  independent  observa¬ 
tions  of  a  decreasing  function.  Maximum  likelihood  estimates  of 
were  therefore  determined  by  using  the  method  of  Brunk  [5]  for  monotone 
functions  applied  to  the  values  of  D^(Zj),  J"1,...,N.  (This  technique 

is  also  used,  e.g.,  in  the  construction  of  maximum  likelihood  estimators 
of  Increasing  failure  rates,  see  [6].) 


Now,  consider  the  following  events: 


(x 

n,n 

<  o,  wx 

>  Cjl. 

Bi" 

iX.  n 

l,n 

a 

o 

A 

>  c2l  * 

*2  ’ 

(X 

n,n 

A 

o 

>  c2J, 

B2- 

[X.  n 
l,n 

>  0,  Wx 

>  C2K 

■ 

n 

< 

(nS  < 

0,  wx  > 

;31» 

B3 

[nS  > 

0,  w  > 

’  n 

;3]. 

J8 


We  know  from  the  symmetry  of  the  normal  distribution  that 


P(\)  =  Pay  =  Dk(ck)  =  e ,  k=l,2,3. 

Hence,  the  procedure  described  above  allows  one  to  determine  the  power 

of  two-sided  tests  and  both  one-sided  tests  on  B. 

Critical  values  for  testing  the  various  hypotheses  of  interest  were 

generated  under  the  assumption  6q  =  B^  =  62  =  8»  The  following  procedure 

2 

is  designed  to  determine  the  power  of  the  tests  for  6q  =  Bj  =  B2  =  £  8. 

First  n  was  fixed  and  N  samples  X^,...,X^  of  standard  normal 
variates  were  generated.  The  following  were  then  calculated  for  each  t 
and  each  j, 


B\  W\  W' 
l  n 


where  the  primed  values  denote  the  fact  that  T^/B  was  replaced  by 
2  2 

T^/ <<5  B).  Since  X^  =  £  (T^/ (6  g)),  one  can  determine  Xj  as  0  function 
of  by  first  solving  for  T^/B  as  a  function  of  and  then 

substituting  this  function  in  the  expression  for  X^,  i*l,...,n. 

Rationaling  the  denominator,  we  obtain 

X|  -  (X1/2)  (6  +  1/6)  +  (6  -  l/6)/l  +  X^/4 ,  i-1 . n. 

Then, 

Ij(k)(6)  -  {A£}  +  {B^},  k-1,2,3, 

were  evaluated  (where  the  prime  indicates  that  primed  values  of  X^, 
i"l,...,n  are  used).  Finally,  the  power  tunctions 
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Pk(6)  =  N  njk)(6)‘  k=1’2*3 

were  computed  as  functions  of  4.  The  values  of  6  considered  were 
6  -  1/5,  1/4,  1/3,  1/2,  1,  2,  3,  4,  5. 

Monte  Carlo  simulations  of  the  power  curves  for  two-sided  tests 
at  e  =  .01  for  k*l,2,3  and  all  values  of  n"2,3,4,5  are  presented 
in  Figures  3,  4  and  5.  It  can  be  observed  that  the  tests  corresponding 
to  k  ■  3  are  more  powerful  than  those  corresponding  to  k  *  1  for 
both  n  *  4  and  n  ■  5  when  testing  Hj  versus  ,  l^  versus  Hj 
or  Hj  versus  H^.  It  is  also  apparent  that  for  fixed  6  the  power 
tends  to  Increase  with  n  ami  that  the  incremental  power  increase  from 
k  ■  1  to  k  ■  3  la  larger  for  n  -  5  than  for  n  -  4. 

On  the  other  hand,  the  aim  ad  power  curves  corresponding  to 
It  “  2 ,  at  c  •  ,01  for  vai  values  of  n,  show  a  markedly  different 

behavior.  One  can  observe  for  two* aided  testH  that  power  first  increases 
as  6  deviates  from  unity  in  either  direction  and  then  decreases  as  4 
approaches  either  aero  or  infinity.  The  tests  are  therefore  not  unbiased 
and  consequently  not  acceptable.  Hence  no  tabulations  of  the  i  (e )  are 
presented.  Tabulations  of  C^c),  C3U)  for  t  -  .005,  .01,  .025,  .05,  .10 
for  n  ■  2(1)5  are  given  in  Table  1  with  the  exception  that  entries  are 
not  possible  for  C^(c)  whenever  c  <  2  n.  Here  *1^(0  is  such  that 
V  5^(0)  ■  c  for  fixed  n  and  k  -  1,3. 
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An  illustration  of  the  use  of  Table  I  is  now  given.  Let 
tj,...,t<.  be  the  ordered  failure  times  of  a  device  in  life  testing. 
Under  the  assumption  that  the  unorderei  .re  times  would  be  a 

random  sample  from  the  two-parameter  distribution  defined  by 
liquation  (0.1),  a  100 ( 1  — t  )%  lower  confidence  bound  on  H  is  given 
by  Fq nation  (3.1),  namely 


min 


5 

\ 

i=l 


5 

\ 

i=i 


2s<r,/v 


A 

Y 


2^3(f)tnS 


A 

V 


Suppose  the  life  times,  in  hours,  were 


t1  =  A8.310 

t2  =  55.15A  t3  =  61,273 

t4  =  58,110  t5  =  67,769 

and  we  set  e  =  .1.  We  compute  a  lower  confidence  bound  for  8  at 
level  .90  to  be  equal  to  A5,375. 
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Figure  3.  Graphs  of  power  as  a  function  of  <5. 
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r 


Figure  4.  Graphs  of  power  as  a  function  of  6. 
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Figure  5.  Graphs  of  power  as  a  function  of 
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